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The Fermi Large Area Telescope (LAT) has made the first definitive GeV detections of the binaries LS I +61°303 
and LS 5039 in the first year after its launch in August 2008. These detections were unambiguous because, 
apart from a reduced positional uncertainty, the 7-ray emission in each case was orbitally modulated with the 
corresponding orbital period. The LAT results posed new questions about the nature of these objects, after the 
unexpected observation of an exponential cutoff in the GeV 7-ray spectra of both LS I +61°303 and LS5039, 
at least along part of their orbital motion. We present here the analysis of new data from the LAT, comprising 
2.5 years of observations through which LS I +61°303 continues to provide some surprises. We find a sudden 
increase in flux in March 2009 and a steady decrease in the flux fraction modulation. The LAT now detects 
emission up to 30 GeV, where prior datasets led to upper limits only. At the same time, contemporaneous TeV 
observations either no longer detected the source, or found it -at least in some orbits- close to periastron, far 
from the usual phases in which the source usually appeared at TeV energies. The on-source exposure of LS 5039 
has also drastically increased along the last years, and whilst our analysis shows no new behavior in comparison 
with our earlier report, the higher statistics of the current dataset allows for a deeper investigation of its orbital 
and spectral evolution. 



1. Introduction 

To date there are only five X-ray binaries that 
have been detected at high (HE; 0.1-100 GeV) or 
very high-energies fVHE; >100 GeV): LS I +61°303 
[EE 9; LS5039 lip!, PSR B1259-63 @, Cyg X- 
3 [J, and Cyg X-l Q. Of these sources, only 
LS I +61°303 and LS 5039 exhibit a non-transient be- 
havior at high-energies and consequently, have been 
detected as persistent, albeit variable, sources of 7-ray 
emission. They also share the property of being (to- 
gether with PSR B1259-63, see sources detected 
at both GeV and TeV energies. The other systems 
mentioned have been unambiguously detected only in 
one band, either at GeV or at TeV. 

The early LAT reports of GeV emission from 
LS5039 and LS I +61°303 were based upon 6-9 
months of survey observations [H, Q. Both sources 
were detected at high significance and were unambigu- 
ously identified with the binaries by their flux modula- 
tion at the corresponding orbital periods, 26.4960 for 
LS I +61°303 [iff and 3.90603 days for LS5039 p3]. 
The modulation patterns were consistent with expec- 
tations from inverse Compton scattering models, and 
were anti-correlatcd in pha se with pre-existing TeV 
measurements (e.g., (T3, Il4j). The anti-correlation of 
GeV-TeV fluxes is in fact a generic feature embedded 
in inverse Compton models describing the TeV fluxes, 
where the GeV emission is enhanced (reduced) when 
the highly relativistic electrons seen by the observer 
encounter the seed photons head-on (rear-on); (e.g., 
[TBI [l6j ) . Fermi measurements provided a generic con- 
firmation for these inverse Compton models. 

Both sources presented exponential cutoffs in their 
high-energy spectra, at least along part of the orbit. 



To be precise, an exponential cutoff was a better fit 
to the SED -as compared with a pure power-law- in 
phases surrounding the inferior conjunction (INFC) 
of LS 5039, and when taking into account the aver- 
age spectrum of LS I +61°303 along its whole orbit. 
Statistical limitations with the used amount of data 
prevented the determination or the ruling out of an ex- 
ponential cutoff in any orbital cut of LS I +61°303 or 
in the superior conjunction (SUPC) of LS5039. The 
spectral energy distributions with the exponential cut- 
offs that were reported were reminiscent of the many 
pulsars the LAT has discovered since launch 17j, al- 
though this was far from a proof of an LS 5039 or 
LS I +61°303 pulsar nature. No pulsations has been 
found at GeV energies. 

In this work we present the results of the analy- 
sis of 2.5 years of LAT survey observations of both 
LS I +61°303 and LS 5039. We investigate any long- 
term flux variations of the sources and explore too the 
possible spectral variability for both systems. We ex- 
pect that our analysis will constrain future theoretical 
studies on these interesting objects. 



2. Observations and data reduction 

The Fermi Gamma-ray Space Telescope, launched 
on 2008 June 11, carries onboard the Large Area 
Telescope (LAT). The LAT is an electron-positron 
pair production telescope, featuring solid state sili- 
con trackers and cesium iodide calorimeters, sensitive 
to photons from -20MeV to > 300 GeV pL{|. It has 
a large field of view with ^2.4 sr (at 1 GeV) and an 
effective area of ~8000cm 2 for >1 GeV. 
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The Fermi survey mode operations began on 2008 
August 4. The full dataset used for this analy- 
sis spans 2008 August 4, through 2010 December 
4. The data were reduced and analyzed using the 
Fermi Science Tools v9r20 package. 1 The standard 
onboard filtering, event reconstruction, and classi- 
fication were applied to the data 18]. The high- 
quality "diffuse" event class was used. Throughout 
the analysis the "Pass 6 v3 Diffuse" instrument re- 
sponse functions (IRFs) were applied. Where re- 
quired in the analysis, models for the Galactic diffuse 
emission (glLiem-v02.fit) and isotropic backgrounds 
(isotropic-iem_v02.txt) were used. 2 

The binned maximum-likelihood method of gtlike, 
which is part of the ScienceTools, was applied to de- 
termine the intensities and spectral parameters pre- 
sented in this paper. We used all photons with energy 
>100MeV in a circular region of interest (ROI) of 10° 
radius centered at the position of LS I +61°303 and 
LS 5039, respectively. For source modeling, the 1FGL 
catalog [i9j . derived from 11 months of survey data, 
was taken. 



3. LS5039 Results 

LS 5039 is located in a complicated region toward 
the inner Galaxy with high Galactic diffuse emission 
and many gamma-ray sources. In an earlier publi- 
cation about this source Q we derived that a power 
law plus an exponential cutoff describes best the data. 
The photon index was r = 1.9 ± 0.1 s tat ± 0.3 sys t; 
the 100-300GeV flux was (4.9 ± 0.5 sta t ± 1.8 syst ) x 
10 _7 phcm~ 2 s _1 and the cutoff energy was found to 
be 2.1±0.3 st at±l.lsystGeV. 

Firstly, we analyzed the orbitally-averaged data of 
LS 5039. As spectral models for LS 5039, we used a 
power law as well as a power law with an exponen- 
tial cutoff, and compared likelihood obtained to test 
the significance of a spectral cutoff. The likelihood ra- 
tio between the power-law and cutoff power-law cases 
clearly indicates that the power-law assumption is re- 
jected. 

Spectral points at each energy band were obtained 
by dividing the data set into each energy bin and per- 
forming maximum likelihood fits for each of them. Re- 
sulting spectral energy distribution (SED) is plotted 
in Figure [TJ 

Following the H.E.S.S. analysis by [lj] and our pre- 
vious analysis, the whole data set was divided into 



1 See the Fermi Space Science Centre (FSSC) 
website for details of the Science Tools: 
http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/ 

z Descriptions of the models are available from the FSSC: 
http://fermi.gsfc.nasa.gov/ssc 



two orbital intervals: superior conjunction (SUPC; 
<f> < 0.45 and 0.9 < (j>) and inferior conjunction (INFC; 
0.45 < < 0.9). The SUPC and INFC data were an- 
alyzed in the same way as the orbitally averaged data. 
Being consistent with our previous paper, the power- 
law assumption for the SUPC spectrum is rejected. 
Although a single power law was not rejected for INFC 
in our previous analysis using 10-month data, a cutoff 
power law is preferred also for INFC with the present 
data set. The corresponding SED is shown in FigurcQ] 



4. LS I +61°303 Results 

LS I +61° 303 is one of the brightest sources in the 
7-ray sky and towers above other emitters in its neigh- 
borhood. This has allowed us to clearly detect, in 
March 2009, a ~30% increase on the average flux. We 
graphically show the flux change in Figure [U by plot- 
ting the folded lightcurves before and after the found 
transition in March 2009. Before the transition, the 
modulation was clearly seen and is compatible with 
the already published phasogram, whereas afterwards, 
the modulation gets fainter. Note that the datasets 
corresponding to the reported results [|| and what we 
here referred to as before the flux change span almost 
exactly the same time range, with the consequence of 
our current analysis essentially reproducing the one 
previously published (3[. The spectra derived before 
and after this flux change are also shown in Figure [2J 
where the increase in flux is obvious. 



4.1. The multi-wavelength context 

LS I +61°303 has also been monitored with the 
i?XTi5-Proportional Counter Array (PCA) and folded 
lightcurves were produced. In the left panel of Figure 
[3l we show a direct comparison between the phaso- 
grams in X-ray and in 7-rays, with simultaneously 
taken data. We divide the whole LAT dataset into 
4 periods of 6 months each and compare them with 
correspondingly obtained PCA X-ray data. The di- 
vision in periods of 6-months is justified in order to 
have enough statistics in 7-rays for each individual 
time bin, and such that orbit-to-orbit X-ray variabil- 
ity does not dominate the flux fraction changes. It 
is apparent that the X-ray modulation is always vis- 
ible in each of these five panels, albeit with variable 
amplitude of flux modulation. Instead, at GeV ener- 
gies, the LAT data indicates that the modulation fades 
away until the variability along the orbit is barely vis- 
ible in the last 6 months of our data. In Figure 2] we 
show the average 7- and X-ray fluxes fitted in each of 
the periods considered. The X-ray flux is increasing, 
whereas the 7-ray flux seems to decrease again after 
the flux change in March 2009. 
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Figure 1: Left: The overall spectrum of LS 5039 over 2.5 years of data are shown in red. Also, very high-energy data 
points from H.E.S.S. are plotted, but one has to keep in mind that they are not simultaneous. The gray filled points 
correspond to inferior conjunction of the system and the open ones to superior conjunction. Right: Spectra for inferior 
(red) and superior (blue) conjunction of the system. Note that all plots are preliminary. 




Figure 2: Left: Folded lightcurve of LS I +61°303 before the flux change (blue), when the modulation is still clearly 
visible. Middle: Folded lightcurve after the flux change in March 2009 (red) when the modulation becomes less 
pronounced. The modulation gets fainter. Right: Comparison of the spectra derived before (blue) and after (red) the 
flux change in March 2009. Note that all plots are preliminary. 



5. Summary 

After two years of continuous data taking of the 
two binary systems LS 5039 and LS I +61°303 we con- 
firm that their energy spectra are best described by 
a power law function with an exponential cutoff. For 
LS 5039 , in comparison with earlier publications of 
this source, with larger data set at hand we can ex- 
tend the observed spectrum to higher energies. No 
significant emission of LS I +61°303 above 30GeV 
has been detected with the Fermi-L AT yet. Further- 
more, in the case of LS I +61°303 variable modulation 



in the folded light curves was observed. Before a flux 
change in March 2009, the modulation was clearly vis- 
ible, whereas afterward it faints away. This behavior 
was not predicted by any theoretical model so far and 
has to be investigated further. 
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Figure 3: Comparison between the 7-ray (blue) and the 
X-ray (gray) data of LS I +61°303 . For each of the 4 
separate 6-months periods, the 3-10 keV (gray) and the 
100 MeV-300 GeV (blue) folded lightcurves are shown. 
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